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Abstract. Integrated ground-based and HST optical 
studies of isolated neutron stars have provided important 
independent datasets in the determination of emission ac- 
tivity, particularly in the fitting of anticipated Rayleigh- 
Jeans extrapolations from EUV/soft X-ray datasets, de- 
spite their intrinsic faintncss. Differentiation of the pulsed 
and unpulsed fluxes and consequently of the nonther- 
mal and thermal modes of emission could provide defini- 
tive data with which to constrain this blackbody contin- 
uum. Based upon high speed photometric observations of 
Geminga and PSR B0656+14 in the B band, we have 
combined upper limits of unpulsed emission with recently 
published model-fits with a view to assessing possible im- 
plications for the R/d parameter. For Geminga, with a 
known distance of ~ 160 pc, we find that Roc < 9.5 km 
for a ~ blackbody source, and Roo < 10.0 km with the 
presence of a magnetized H atmosphere. In addition, we 
suggest that PSR B0656+14 is some ^4-5 times closer 
than the 760 pc estimated from DM measurements alone. 
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1. Introduction 

Uncertainty regarding the behaviour of nuclear matter in 
the deep neutron star interior has compromised a complete 
description of the dense matter equation of state (EOS). 
Various theoretical models in circulation predict a range 
of macroscopic observables, such as masses, radii, temper- 
atures and maximum rotation rates, and as such are open 
to scrutiny with empirical data. Relativistic effects com- 
plicate such observations, with the true radius, R, related 
to the apparent radius at a distance ( d — > oo ) as 



R r , 



(1 + z)R = R/y/1 - 2GM/Rc 2 , 



(1) 



with M being the mass of the compact object, and z the 
associated gravitational redshift. Theory to date indicates 
that 7km < R^ < 20km, dependent on the stiffness of the 



EOS (e.g. Lindblom 1992). Binary studies suggest a com- 
mon value of neutron star mass approximately that of the 



canonical estimate of 1.4 M Q (Van Kerwijk et al. 1995). 
Typically, the model neutron star is assumed to have M 
~ 1.4 M Q and R x ~ 13km (R - 10 km). Based on ob- 
servations of the neutron star's thermal emission in the 
extreme UV (EUV) and soft X-ray bands, it should in 
principle be feasible to compute the ideal blackbody spec- 
tral energy distribution (SED) as a function of (T sur f ace , 
TVh & R/d), and for a known distance d, an estimate of 
the apparent neutron star radius. It is generally agreed 
that such a measurement would have a profound impact 
in constraining EOS models. However, strong galactic HI 
absorption at these wavelengths restricts observations to 
the closest neutron stars, and uncertainty in X-ray detec- 
tor sensitivities at these low (< 0.2 keV) energies ( e.g. 
Walters & An 1998) has compromised attempts to accu- 
rately determine estimates for Ah and T sur f acc . Further- 
more, such an analysis may be complicated by neutron 
star phenomenology, such as atmospheric opacity effects, 
an active magnetosphere, hot polar cap regions and ac- 
cretion processes. It is critical to separate the various 
contributions, so as to estimate the genuine total surface 
component. For the older, isolated neutron stars (INS), 
phase-resolved studies in the soft X-ray, EUV and opti- 
cal wavebands can provide such a SED, and thus a real 
possibility of determining R^ . The Rayleigh- Jeans tail in 
the optical regime provides stringent constraints to any 
SED model- fit, and in this waveband atmospheric opac- 
ity effects are expected to have the most noticeable im- 
pact. In the X-ray regime where the blackbody contin- 
uum peaks, low Z atmospheres preferentially transmit 
radiation from the lower, hotter regions of the photo- 
sphere, producing a X-ray spectrum suggestive of a 'hot- 
ter' source ( Romani 1987 ) . The Rayleigh- Jeans tail is un- 
affected by this deviation, and discrepancies between op- 
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tical spectrophotometry and higher energy extrapolations 
may be rigorously tested (Pavlov et al. 1998). However 



the intrinsic faintness of these astrophysical objects in 
the optical regime coupled with the limitations of cur- 
rent technology have restricted previous optical stud- 
ies to deep integrated photometry. These observations 
have in some ways aided such ther mal continuum stud- 
ies (e.g. Walters fc Matthews 1997 ), but the differentia- 
tion of pulsed (predominantly nonthermal) and unpulsed 
(thermal) components would be ideal in a more rigorous 
treatment. Recently, the TRIFFID high speed optical pho- 
tometer detected pulsations in the B band from the optical 
counterpart s of the middle aged pulsars Geminga an d PSR 
B0656+14 ( |Shearer et al. 19971 [Shearer et al. 1998j ). Both 
light curves are highly pulsed and suggest a dominant non- 
thermal mode of emission optically. Despite their extreme 
faintness, it was possible in both cases to determine up- 
per limits to each pulsar's thermal component of emission. 
In this letter, we combine these unpulsed limits with the 
results of recently published ground-based and HST pho- 
tometric analysis on these two pulsars ( [Pavlov et al. 1997 , 
Martin et al."T99g| hereafter PWC97 k MHS98), and de- 
rive radius/distance estimates for both. In particular, for 
the case of Geminga with a known parallax distance of ~ 
160 pc, we provide for the first time an upper limit for 
Roc based upon such phase-resolved photometry. 

We approach the analysis of the optical photometry 
by assuming a model fit incorporating both nonthermal 
power-law and thermal components of emission in the UB- 
VRI regime, as originally adopted by PWC97. This two- 
component model fit is defined as: 



/(") 



/o(- 



v 

9o 1 — 



(2) 



with vo = 8.766 x 10 Hz an arbitrary reference frequency, 
and A(y) the interstellar extinction determined following 
Savage & Mathis (1979). go and fo are taken to be the 
values of the nonthermal and thermal fluxes at v = z/q. 
The latter can be expressed, for the chosen value of vq, as 



go = 3.116 x 10~ 31 G 



erg 



cm 2 s Hz 



G = T 6 



R 



in 



-tson 



(3) 



where T = 10 6 Tg K is the apparent neutron star bright- 
ness temperature, Rao — 10i?io km the radius and d = 
500d5on pc the distance to the blackbody. Model fits for 
both pulsars have suggested the presence of both emission 
modes, yet uncertainty in the optimum (T sul f aco ,iVH, R/d) 
solution based on EUV/soft X-ray datasets restricts accu- 
rate differentiation. By providing upper bounds for the 
unpulsed flux, and indirectly N-r via the estimated A(y) 
towards the pulsar, we can derive an independent estimate 
of G and in this way, constrain (T sur f ace , R/d) space. The 
use of the lower bounds to T sur f ace for a given neutron star 
would then yield upper bounds to the R/d parameter. 



2. Geminga 

PSR 0633+17, or Geminga, provides one of the most ideal 
neutron stars for such an analysis, as HST observations 
have determined its distance by parallax to 1 59+34 pc 
flCaraveo et al. 1996| ). However, it is by no means clear 
precisely what the nature of its emission processes are, 
despite considerable observational efforts from the op- 
tical to 7-rays. Original ROSAT observations suggested 
dominant thermal emission from the surface modulated 



by a hotter polar cap (Halpern & Ruderman 1993), al- 



though later observations especially those of ASCA indi- 
cated the emission was thermal and nonthermal in ori- 



gin ( Halpern fc Wang 1997 ). Disagreement remains on 
the optimum SED fits to data from ROSAT, EUVE 
and ASCA with the pulsar's surface temperature within 
(2.5 — 6) x 10 5 K. There have been further sugges- 
tions that the pulsar's thermal emission is strongly af- 
f ected by magneto spheric cyclotron resonance blanketing 
( Wang et al. 1998 ) , which would severely compromise any 
estimate of the true neutron star surface temperature. Op- 
tical observations suggest an unusual SED, with devia- 
tions in the expected functional form, possibly as a re- 
sult of ion cyclotron absorption/emission processes. The 
most recent of such observations (MHS98), made using 
the Keck LRIS, spanned 370-800 nm and yielded a flat 



power law shape (f v oc 



.-0.8 



) and a noticeable broad dip 



at 630-650 nm with perhaps a slight modulation at V, B 
& / as advanced by Bignami et al. (1996). The compos- 
ite power-law may be fitted by either a combined black- 
body and power-law (i.e. nonthermal emission with a ~ 
1.9 ± 0.6) or a blackbody plus global ion cyclotron emis- 
sion. The former model is undoubtedly the most likely, 
following the discovery of a highly pulsed lightcurve in the 
B band (Shearer et al. 199S). Fig. 1. shows the model-fit 
spectrum of MHS98 with the results of Shearer et al. 
(1998). 

A value of G — 4.4 is obtained by applying a so- 
lution in the form of (3) to the best-fit Keck obser- 
vations of MHS98, setting A(i/) ~ due to the pul- 
sar's close proximity. Introducing our independently es- 
timated unpulsed upper limit in the B band of 6.77 x 
10~ 31 ergs cm 2 s _1 Hz -1 , we conclude with G < 3.6±0.1, 
allowing errors in G for the spectral response of the 
MAMA/S-filter combination. The decrease can be in- 
ferred as either a reduction in the emission area or a drop 
in T sur f aco (d ~ 160 pc). As a consequence, one expects 
the nonthermal component to increase proportionately in 
the original model fit, if this unpulsed estimate is rep- 
resentative of a general surface based thermal emission. 
If we assume that ^surface remains unchanged, then the 
obvious conclusion is that R has been overestimated. In 
fact, using G < 3.6, T surfaco - 4.0 x 10 5 K (MHS98) and 
d = 159jl|| pc suggests R^ < 9. 5~tf,o km for a — black- 
bo dy source (indisting uishable from a Fe/Si atmosphere - 
see Pavlov et al. 1998). An alternative explanation for the 
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Fig. 1. Keck LRIS Optical spectrum of Geminga, with best- fit 
two component model of Martin et al. (1998) and TRIF- 
FID/BTA pulsed and unpulsed flux estimates in B (Shearer 
et al. 1998) 



change in G could be understood under the assumption of 
atmospheric opacity effects for the total surface thermal 
emission at X-ray energies. As PWC97 point out, there 
is no ideal way to reconcile the actual brightness temper- 
ature with the estimated X-ray T sur f acc . Fits with magne- 
tized ( ~ 10 12 G ) H atmospheres suggest that a spectrum 
in the optical regime can be fitted as a Rayleigh- Jeans tail 



rfacc 



with TL, 



yields < 



lO.OtH km. 



3. PSR B0656+14 

An extensive literature exists devoted to the analysis of 
EXOSAT, ROSAT and EUVE observations of the pul- 
sar PSR B0656+14. The consensus is of thermal emis- 
sion from the surface, modulated by emission from a hot 
polar cap and from nonthermal activity in the magne- 
tosphere. However, its distance estimate is by no means 
exact ranging from ~ 100 to 760 pc ( paraveo et al. 1994 ) 
- the latter based upon a DM fit to the uncertain galactic 
electron model in this vicinity (Taylor et al. 1993). Ap- 
plying a common distance estimate of ~ 500 pc, a range 
of solutions are possible in (T sur f acc , Nr) space assum- 
ing the canonical model, and this is reflected in the lit- 
erature, with r atface = (3 - 9) x 10 5 K and A% = (0.5 
2.0)xl0 20 cm 2 (|Finlcy et al. 1992j |Anderson et al. 1993| 



Grciveldingcr et al. 1996). Early limite d observations us- 
i ng the NTT ( I Caraveo et al. 19941) and the HST 



( Pavlov et al. 1996 ) indicated that the optical counter- 
part's emission was predominately nonthermal in nature, 
and a subsequent two-component model fit to detailed 



HST and ground-based photometry spanning the UB- 
VRI regime by PWC97 substantiated this assessment - 
although the detection of highly pulsed emission in the B 



band (Shearer et al. 1997) from the optical counterpart 



unequivocally confirmed this hypothesis. The upper limit 
on the unpulsed flux from the resulting lightcurve was es- 
timated to be 8 x 10~ 31 ergs cm -2 s^ 1 Hz -1 . Fig. 2. shows 
the PWC97 model fit with these pulsed/unpulsed fluxes. 

PWC97 fitted the observed UBVRI spectrum with a 
two-component model following the formalism of (2) and 
(3). Applying the interstellar extinctions determined for 
the three estimated colour excesses towards the pulsar 
of E(B - V) = 0.01, 0.03 and 0.05 yielded best fit val- 
ues of G = 3.0, 3.7 and 4.3 respectively. Taking g — 
8 x 10~ 31 ergs cm -2 s _1 Hz -1 , (2) may be rearranged in 
terms of G, and solved for the three E(B — V) estimates, 
concluding with G = 4.4, 4.8 and 5.2 (all ±0.1). respec- 
tively. An increase of this G parameter is consistent with 
either an increase in the expected Tsurfacc, emission area 
or a decrease in the pulsar's distance. 



0.9T off (PWC97). Applying this correction ~ 



HST & BTA Best-Fit (Pa 



p = 25.1±0.3 
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Fig. 2. HST & 6m BTA based photometry and best-fit two 
component model of PSR B0656+14 for the colour excess 
E(B - V) = 0.03 (Pavlov et al. 1997), with TRIFFID /BTA 
pulsed and unpulsed flux estimates in B superimposed (Shearer 
et al. 1997) 



In terms of previous model fits, G ranges from 0.9 
flFinlcy et al. 1992}), 2.1 (preiveldinger et al. 1996|) to 2.6 



( Anderson et al. 1993[ ) , the latter incorporating a magne- 
tiscd (~ 1Q LZ G) H atmosphere. As observations have 
accrued, and uncertainties in X-ray detector sensitivities 
have been addressed, the trend has been a decrease in the 
derived surface temperature, from ~ 9.0 x 10 5 K to ~ 5.0 x 
10 5 K. Recent work by Edelstein et al. (1998), which sub- 
stituted EUVE DS data in place of the uncertain low en- 
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ergy ROSAT PSPC channels, has yielded a (T surfaco ,Af H ) 
space differing markedly from earlier ROSAT results 
alone. Combining this new parameter space with inde- 
pendent model fits incorporating the B impulsed upper 
limit suggests T surfaco > 5.0xl0 5 K QGolden 1998|). It is 



possible to constrain R/d for B0656+14 in two ways - 
firstly by determining apparent expected radii using both 
column density and DM derived distances, and secondly 
by applying canonical or otherwise derived radius limits 
to yield optimum distance scales. Nn estimates towards 
PSR B0656+14, although by no means certain, suggest 
that the pulsar is y 250-280 pc ( Anderson et al. 1993} , 



Edclstcin ct al. 1998 ), rather closer than the DM derived 
distance of 760 ± 190 pc. By manipulation of (3) with 
^surface > 5.0 xlO 5 K and the range of G parameters de- 
termined from the unpulsed upper limit, the expected ra- 
dial estimates for the JVh distances are 14.7 < Roo < 17.7 
km, and substantially greater in the case of the radio de- 
rived distance. We note that such radii estimates are in 
excess of the 14 km upper limit determined by Walter et 
al. (1997) for the old INS RXJ185635-3754. Alternatively, 
applying the ideal canonical Roo ~ 13km, manipulation of 
(3) as above suggests that 205 < d < 227 pc based upon 
the range of E(B — V) estimates. This supports the con- 
clusions of PWC97 although placing the pulsar somewhat 
in closer proximity than had been originally thought. In- 
deed if one was to consider the proposed estimate of Roo 
~ 9. 5^2 o k m f° r Geminga as a working upper limit, this 
would place PSR B0656+14 at a distance of no less than 
d = 152^32 for the suggested optimum colour index of 
E{B - V) = 0.03 (PWC97) 



4. Conclusions 

The acquisition and analysis of optical photometric data 
on INS has been shown to provide an independent dataset 
from which constraints to the optimum thermal spectral 
energy distribution may be applied. We have attempted 
for the first time to apply the phase-resolved optical flux 
consistent with that expected specifically from the ther- 
mal component to the previous integrated analysis of the 
pulsars Geminga and PSR B0656+14. Despite being up- 
per limits to the unpulsed optical flux in the B band, we 
find that in both cases, the resulting blackbody spectral 
distribution is constrained to the extent that we can set 
upper limits to the R/d parameter for both neutron stars. 
For the case of Geminga, with a known parallax derived 
distance of ~ 160 pc and using the lower Tsurfaco limit, 
we suggest that Roo < 9.5 km for a ~ blackbody source, 
and Roo < 10.0 km with the presence of a magnetized 
H atmosphere. Previous work using these unpulsed upper 
limits has suggested a T sur f aco > 5.0xl0 5 K for the pulsar 
B0656+14 (Golden, 1998), and under the assumption of 
Roo ~ 13km, places the pulsar at ~ 210 pc - in contrast to 
the DM derived distance of 760 pc. Assuming the neutron 
star has Roo < 9.5 km, then this limits the distance to ~ 



160 pc. This suggests the possibility that the pulsar may 
be a viable candidate for a parallax measurement attempt 
using the HST. 

Despite using upper limits to the unpulsed fluxes of 
these two INS, we have been successful in indicating how 
important such definitive measurements are in the rigor- 
ous derivation of a given neutron star's thermal parame- 
ters. In this way, we have independently provided limits to 
the radius of one, Geminga, and the distance of another, 
PSR B0656+14 and shown the promise to future studies 
this observational technique offers. This short report doc- 
umenting these results is most timely, following on from 
a recent forum on the study of INS in this energy regime, 
at which the future importance of such high speed pho- 
tometric studies was stressed (Romani 199S), and where 
a recent theoretical analysis based upon the discovery of 
RXJ185635-3754 has shown that neutron star radii should 
be ideally < 10 km if one can hope to effectively constrain 
the EOS (An et al. 1998). That our results point to such 



a radius limit only emphasises the need for further high 
speed photometry of both pulsars, particularly Geminga, 
so as to provide definitive, rather than upper limit, fluxes 
of their unpulsed emission. 
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